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Recently there has been a renewed interest in the chemical physics literature of factorization of the position 
representation eigenfunctions {$} of the molecular Schrodinger equation as originally proposed by Hunter in 
the 1970s. The idea is to represent $ in the form ipx where x is purely a function of the nuclear coordinates, 
while tp must depend on both electron and nuclear position variables in the problem. This is a generalization of 
the approximate factorization originally proposed by Born and Oppenheimer, the hope being that an ‘exact’ 
representation of <i> can be achieved in this form with ip and x interpretable as ‘electronic’ and ‘nuclear’ 
wavefunctions respectively. We offer a mathematical analysis of these proposals that identifies ambiguities 
stemming mainly from the singularities in the Coulomb potential energy. 


I. INTRODUCTION 

In the Born-Oppenheimer (BO) model an eigenfunc¬ 
tion of the molecular Hamiltonian is often approxi¬ 
mated by a simple product of ‘electronic’ and ‘nuclear’ 
wavefunctions 

$(r, R) Ki p{r, R)xiR), {T\T)r = 1 for all R (1) 


with the nuclear function /(i?) defined as a marginal by 
means of 

|/(R)|2 = y$(r,i?)>(r,i?)dr. (3) 

This fixes f{R) to within a phase factor, 

(4) 


where p is an eigenfunction of the clamped-nuclei elec¬ 
tronic Hamiltonian, and x is a vibration-rotation wave- 
function for the nuclear motion on the potential energy 
surface (PES) associated with p. r and R stand collec¬ 
tively for the electronic and nuclear coordinates respec¬ 
tively. This is probably the most important approxima¬ 
tion in the quantum theory of molecules, and it has been 
subject to exhaustive mathematical analysis^. It should 
be noted that the BO approximation is a singular per¬ 
turbation problem based on the limit of inhnite nuclear 
massi^— , and it is plausible that the product form, when 
valid, is a direct consequence of the limit. Nevertheless 
there has recently been renewed interest^^— in the ques¬ 
tion as to whether the ~ symbol in m can be replaced 
by = with suitably redefined ‘electronic’ and ‘nuclear’ 
functions p and x, a so-called ‘exact’ factorization of an 
eigenfunction $. The present work offers a mathematical 
analysis of this factorization. 

The earliest attempt that we know of to write an exact 
wavefunction $(r, R) of the Coulomb Hamiltonian H for a 
neutral system of electrons and nuclei in a factorized form 
was made by Hunter—. His analysis was made in terms 
of conditional and marginal probability amplitudes; the 
wavefunction $, assumed normalized, is written in the 
product form 


<L>{r,R) = f{R)cj){r,R) (2) 


Manuscript published October 13, 2015: J. Phys. A: Math. 
Theor 48 (2015) 445201 (20 pages) 


where 0{R) is a real function of the nuclear coordinates. 
In the absence of a criterion to choose it. Hunter sug¬ 
gested taking simply 


f{R) = -\fiR)\or f{R) = \f{R)\. (5) 


The associated function, (j), is then defined to be the 
quotient 


cj}{r, R) 


$(r, R) 

~1W 


( 6 ) 


and it satisfies the normalization condition 
J (j){r, R)*(j){r, R)dr = 1 

for all R. Obviously this construction would be problem¬ 
atic if either f{R) has zeroes for hnite R or $(r, R)/f{R) 
is too irregular at infinity. The construction, if appli¬ 
cable, is available for any wavefunction and so the 
nuclear functions {f{R)i} are required generally to be 
quite different from the usual approximate nuclear wave- 
functions for vibrationally excited states which do have 
nodesiS. Equally, it is evident that every wavefunction <I>i 
has its own distinct electronic factor, (j)i; this is to be con¬ 
trasted with the BO description dl} where whole groups 
of approximate eigenfunctions {$„} share one electronic 
state p which supports a vibration-rotation manifold of 
states. 

More recently, attempts at factorization have been 
made by Gross and co-workers^“— using both time- 
independent and time-dependent formulations which are 
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a development of Hunter’s approach and, reverting to 
the time-independent form, Cederbaum has proposed a 
related factorization scheme^ (see also the subsequent 
Erratum^); it is these more recent proposals we anal¬ 
yse here though we point out the connection between 
Hunter’s early work^ and the more recent work of Gross 
et ai. 

For a freely moving system it is always possible to sep¬ 
arate completely the centre-of-mass dynamics (free mo¬ 
tion) from the internal motions of the molecule. It is not 
essential for an account of factorization to remove the 
centre-of-mass motion (and to do so complicates the form 
of the internal Hamiltonian somewhat - see below) but if 
one does not, the description of the bound-states of the 
molecule is more involved. If we denote the position coor¬ 
dinate of the centre-of-mass by R, and introduce a set of 
independent internal position variables, straightforward 
calculatioiiii yields the Hamiltonian H in Schrddinger 
representation separated into internal and centre-of-mass 
contributions 

H = , H' = H,z+T„. (7) 

Here Mt is the total molecular mass and Hez, which 
only differentiates in the electronic variables, and Tn 
are defined explicitly in Section [TTl The spectrum of the 
Hamiltonian H is purely continuous and the description 
of the molecular bound-states requires a rather compli¬ 
cated mathematical formulation. Indeed, H has no eigen¬ 
functions. On the other hand the eigenfunctions of the 
internal molecular Hamiltonian H', (1171) . are square inte- 
grable and have (distributional) derivatives up to second- 
order that are also square integrable; they belong to the 
Sobolev space These are the true bound-states {^'i} 
associated with discrete energies {Ei} that describe the 
internal motions of the molecule. The regularity prop¬ 
erties of the bound-state wavefunctions turn out to be 
important for a precise account of factorization, and it is 
advantageous to frame the discussion directly in terms of 
the internal Hamiltonian. 

The paper is organized as follows; in Section |IT] we re¬ 
view some features of the Schrddinger equation for the 
Coulomb Hamiltonian which are pertinent here. We em¬ 
phasize the occurrence of singularities in the Coulomb po¬ 
tential energy operator which require that the molecular 
Schrddinger equation (in position representation) be in¬ 
terpreted in a more general setting than a classical partial 
differential equation. Factorization of an eigenfunction of 
the Coulomb Hamiltonian evidently must take account 
of this mathematical setting, but it also brings in new 
problems which this paper aims to characterize. In the 
following we shall identify some ambiguities in the recent 
discussions of factorization and describe a mathematical 
framework where the ambiguities are removed. Some re¬ 
maining uncertainties are reported. A common feature 
of the proposed factorization schemes is that the factors 
(^, x) are solutions of a system of non-linear equations. 
In Section lHll we study the formal computation that gives 


rise to these equations, paying close attention to the reg¬ 
ularity properties required to give it a precise meaning. 
In Section IIVI we discuss a variational calculation that is 
related to the system of non-linear equations. 

Two particular facts are important for any proposed 
factorization. Firstly, even if one can find factors (^, x) 
both of which belong to the appropriate TL^ Sobolev 
spaces, in the r and R variables respectively, it does not 
follow automatically that their product m belongs to 
TL^ (in all variables r,R). Secondly it is possible that 
an electronic function (/)(r, R) defined in the manner of 
Hunter, dH]), does not belong to the Sobolev space 
and so cannot be interpreted as a bound-state electronic 
wavefunction. We give a model example in Appendix 
IVHIDI Thus both approaches involve technical difficul¬ 
ties that must be overcome and both require a supple¬ 
mentary check that the product Tpx is actually an eigen¬ 
function of the Schrddinger equation. 

An alternative to Hunter’s interpretation of the quo¬ 
tient 4){r, i?), which we will also explore here, is to regard 
it as a molecular wavefunction required to belong to the 
Sobolev space 'H'^{r,R). An example of such a factor¬ 
ization which largely avoids the troublesome technical 
details is given in Section |Vl Although it lacks the di¬ 
rect physical interpretation of Hunter’s approach it does 
incorporate precise information about the behaviour of 
eigenfunctions at infinity. In Section IVII we discuss the 
normalization of the factors proposed in the light 

of the previous results, and correct the modified com¬ 
putation in^; the relationship between the ‘electronic’ 
and ‘molecular’ wavefunction interpretations of <(>, equa- 
tion(l6]), is discussed here. Finally we try to draw together 
our findings in Section IVHI An Appendix reviews some 
key mathematical notions in an informal way. 


II. THE COULOMB HAMILTONIAN 


A molecule considered as a quantum mechanical col¬ 
lection of electrons and nuclei is customarily described 
by the usual Coulomb Hamiltonian H acting on an Eu¬ 
clidean configuration space defined by the particle coor¬ 
dinates with Schrddinger equation 

H$ = E^. (8) 


More explicitly, for a system of W electrons and iV„ 
atomic nuclei we have 


H = 



(9) 


The configuration space is K3Af„-i-3Afe^ where R denotes 
the set of real numbers. The interparticle distances are: 
r-- = |x- - x-|,rf; = |xf - x"|, and rf/ = |xf - x-| 
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in terms of particle coordinates {x®,Xg} in a laboratory 
frame. Here and elsewhere we use (*,j) and {g,h) as 
indices for electrons and nuclei respectively. 

The primes on the summation symbols mean that 
terms with identical indices (‘self-interactions’) are to be 
omitted. We define the set of nuclear collisions as the 
set of those configurations such that = 0 for some 
g,h. Similarly, the set S of all collisions is the set of those 
configurations such that = 0 for some g,h, or r™ = 0 
for some i, h, or = 0 for some j. Collisions have im¬ 
portant consequences for the analytical properties of the 
eigenfunctions which seem to have been first considered 
by Katoi^. In particular, it is expected that cusps appear 
in the wavefunction at collisions if the wavefunction does 
not vanish; it is also possible for an exact eigenfunction 
to have a node at the singularit} d^d"^ . Whether the inclu¬ 
sion of such cusps by means of trial wavefunctions involv¬ 
ing rf? in electronic structure calculations could improve 
their accuracy has been quite widely studied; an example 
can be found ir>i^. We shall consider the matter further 
from a mathematical standpoint!^ when examining the 
proposed factorisations. 

For discussions of molecules where one wants to focus 
on the approximate separability of electronic and nuclear 
motions, it proves convenient to make a specific choice of 
the internal coordinates. The nuclear position variables 
im can be chosen as a set of Nn — 1 translationally in¬ 
variant variables, defined in terms of the original nuclear 
position coordinates, such that one of the new variables 
is the position coordinate for the nuclear centre-of-mass 
X. The electronic coordinates {r®} are a set of W vari¬ 
ables defined in terms of the original electronic coordi¬ 
nates byii 


X. = r. 


X. 


( 10 ) 


is 


^2Nr.-l 

H"(R")=_ ^ _V(R^).V(R)() 

^ g,h=l 

^ ^ ' Z,z, 


(14) 


where rg^(R") is the internuclear separation distance ex¬ 
pressed in terms of the {R”} and the inverse mass matrix 
1 /is in standard formii. 

The electronic and nuclear motions are coupled only 
via a potential term, 


Hen(Rn,re) 



Afn We 

EE 


5=1 


j,en 

39 


^ _ 

(r%R") 


(15) 


where the electron-nucleus distance expression |x® — 
Xg I = is again expressed in terms of the internal co¬ 
ordinates. In the following it will be convenient to write 
{r®} as r, and {R"} as R for simplicity, and denote the 
gradient operator on nuclear coordinates as V”. The first 
(sum) term in (1141) is the kinetic energy operator T„ in 
0; we write it in this shorthand notation as 

T„ = Ev".V'*. (16) 

zg, 

The operator Hg; in ([7]) is composed from the sum of the 
terms in (ITTl) . (ITSl) and the last term in (fHl) . 

The Schrodinger equation for H' defined by (fTTI) - (fT5l) 
is formally an elliptic partial differential equation (PDF) 
in the coordinates (r, i?), 

H'T = E'S, (17) 


With this choice of coordinates the translationally in¬ 
variant Coulomb Hamiltonian takes the form, 

H' ^ H^(r^) + H'*(R") -f H^’"(R",r®). (11) 

The part of the Hamiltonian which can be associated 
with electronic motion is 


^2 We 


2M]sf 




We 


E 


Sttco |r^ - rj 

i,j=i J 


( 12 ) 


where 


Nn 

Mm = 

3=1 


(13) 


The part that can be associated with nuclear motion 


on the reduced configuration space X = 

The occurrence of the Coulomb singularities in H' and 
the physical interpretation of ik require that dni must be 
placed in a more general mathematical setting involving 
the notion of distributional derivatives if it is to be given 
a precise meaning; we refer to Appendix IVIII Al for details 
on distributional derivatives. One has to view (ITTI) in the 
following way. 

Let us denote by C?{X') the set of square integrable 
functions on X. We define the Sobolev space 'H^(A’) 
as the space of £^(A’)-functions such that their distri¬ 
butional derivatives up to second order all belong to 
C3{X'). For T G 'H‘^{X), each term in (fT71) makes sense 
as a £^(A’)-function and the equality takes place in this 
space C?{X). For instance, the term T^ik is a C'^{X)- 
function that satisfies, for all smooth functions h on X 
with bounded support. 


(Jn'^\h) c'^(X) — ( 18 ) 

where T„/i is now computed in the usual sense. One 
can see an eigenfunction ik G L‘^{X) as a distributional 
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solutionto^. This means that for all smooth functions III. EIGENVALUE EQUATION VERSUS NON-LINEAR 
h on X with bounded support, SYSTEM 


(19) 

Essentially what is done here is the differentiations in 
H' are transferred to suitably smooth functions h, us¬ 
ing integration by parts, as required. Note this point of 
view is already necessary in the simplest case: the Hydro¬ 
gen atom. After removal of the centre-of-mass motion, 
the internal Coulomb Hamiltonian involves the electron- 
proton relative coordinate r. The groundstate is given 
by tpo = cexp(—|r|), in appropriate units. This func¬ 
tion is continuous everywhere, and differentiable outside 
the collision at 0. But it is not differentiable at 0 and 
so C21) cannot be understood in the usual way. Now if 
the potential energy terms were smooth functions, for 
example Hooke’s Law for coupled oscillators, the refor¬ 
mulation just described would yield (smooth) solutions 
everywhere that were solutions of the PDE (flTl) in the 
usual sense. It is the occurrence of the singularities in the 
Coulomb potential that cause the main difficulties (to be 
discussed below) for the idea of an ‘exact factorization’ 
of a molecular wavefunction. Eor future reference we de¬ 
note 

as the Sobolev space for wavefunctions 
depending on only the Nn — 1 nuclear coordinates; it is 
contained in the corresponding space 

The best known regularity of an eigenfunction dt of 
the Coulomb Hamiltonian H' is only that its first distri¬ 
butional derivatives are bounde d^^’^^ ; in particular, we 
do not know if one can differentiate everywhere in the 
usual sense. Fortunately, we have some further informa¬ 
tion on d'; we know by elliptic regularityii (see Appendix 
IVHI Gl) . that the following two statements are valid: 

tk is a real analytic function outside the set E. (20) 
tk has at most isolated zeroes outside E. (21) 


In this Section we study the factorisation of eigenfunc¬ 
tions and the associated non-linear problem solved by the 
factors that were presented by Cederbauro^i^ and in the 
contributions of Gross et ah^. We also review Hunter’s 
factorisation^. To begin with we follow Cederbaum’s ar¬ 
guments applied to an eigenfunction of H'. 

Firstly recall that Hunter started not from an assump¬ 
tion of ‘nuclear’ and ‘electronic’ factors but from an ex¬ 
act wavefunction for the molecular system, which he then 
analysed in terms of conditional and marginal probabil¬ 
ity amplitudes to yield a factorization. Cederbaum’s ap¬ 
proach is rather different. In his equation (7a) it is as¬ 
sumed^ that a product form ip^r, R)x{R), where these 
functions are the putative solutions of a pair of coupled 
equations, can represent an exact wavefunction 4', rather 
than that the exact wavefunction can be written in prod¬ 
uct form; equation (7a) thus needs an existence proof. 

Starting from a normalized solution 4^ of H'4> = 
A4' and making the ansatz that 4' (r, R) factorises into 
Tp(r,R)x{R)^ one can try to follow Cederbaum’s discus¬ 
sion (§HAii), disregarding for the moment the question 
of normalization that will be studied later in Section lyll 
Cederbaum’s formal computationi^ can be summarized 
as follows: 

0 = {H'-E)Tpx 

0 = ^(T„x)--V"x-V"^ + x(H'-A)^ (22) 

0 = x(x-'(TnX)^--x“'V”x-V"^ 

+ (H' - E)^) (23) 

° = (H'-^X-'V"x-V"-Ae,(i?))^ (24) 

with 


If we replace the Coulomb interaction by some smooth 
potential, then elliptic regularity shows that 4' is smooth 
everywhere. This explains why the derivatives in (1171) 
can be computed in the usual sense in such a case. 

These facts about a Coulomb eigenfunction already 
yield some useful information about its putative factors. 
Thus, for example, if 4' can be written as x(.R )^(''’5 .R)) 
then neither x nor ||^||r := (^, ^)r can vanish outside E„. 
Indeed, if x or \\Tp\\1 vanishes at some i?o then so does 
4' on {(r, A); i? = i?o}. If i?o ^ then this set con¬ 
tains at least a segment outside E. This contradicts the 
fact, m, that 4' has isolated zeroes outside E. On the 
other hand since we do not know if di* has usual deriva¬ 
tives everywhere we cannot reasonably assume that the 
factors Tp and x nre everywhere regular. This informa¬ 
tion is directly relevant to our consideration of a system 
of equations that provide a formal definition of factors Tp 
and X- 


Eel{R) =-x-\TnX)+E. 
We can rewrite (HI as 


Help} = Eel{R)(p 


with 


Hei = H'-^V"x-V" 


(25) 


(26) 


(27) 


Multiplying (E51) by tp* and integrating over the electronic 
variables, we get 


Eel{R)\\p>\\l = {ip\Helip)r- 
From (1^ and (1^ . we derive 


TnX 





(28) 


(29) 
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Now, if (^, x) solves the coupled, non-linear equations 
(l26)l and (l29)l . then, reversing the above computation, 
we get (|^ and ^ is a solution of (H' — E)Tpx = 0. 

The present paper does not offer a detailed investiga¬ 
tion of the non-linear equations (l26l) and (l 2 ^ but stud¬ 
ies their relationship with the factorization. What is re¬ 
quired is a framework in which the above computation 
can actually be realized. A natural assumption to make 
would be that Tp g 'H^(A’) and x G since 

one wants to interpret them as wavefunctions. Further¬ 
more their product is to be an eigenfunction, Ipx = 
and di* g is therefore essential. Since the lat¬ 

ter property is not guaranteed by ^ g 'H^{X) and 
X G (see Appendix IVIII BIl one could try 

to study the non-linear equations in a subspace "Ho of 

n^{x) 

no = mx)&n\x)xn\M.^^--^y,ipx&n\x)}. 

(30) 

On the other hand it is not at all obvious that all of 
the steps ((Ml) - (l29l) in the computation are valid in the 
setting O- 

First of all, since x may have zeroes, the meaning of the 
division by x hi the above formulae requires explanation. 
This is actually a delicate issue; to see this, let us take an 
example. Let / be a smooth function on X with support 
in the region {|r| < 1; |i?| < 1} and such that / = 1 for 
(r, A) close to 0. Let 5 be a smooth function on 
with bounded support such that g{R) = 1 for |i?| < 2 . 
Consider the smooth function x{R) = 5 (i?)|i?p™, for 
some integer m. Away from i? = 0, //x is a smooth 
function (given by /(r, but it is a quite nasty 

function near i? = 0 if m is large enougbi^. In particular, 
it does not belong to C^{X) and it is not clear how to 
interpret it as a distribution. 

Thus we must consider how we might give a definite 
meaning to the computation (IMl) - (l29l) . In our first ap¬ 
proach which is a local treatment, we avoid the set E, 
that is we restrict the values {r, i?} to lie outwith S. In 
view of (IMl) it would be natural to assume that the fac¬ 
tors Tp and X have at least conventional derivatives up to 
second order in this region. Since both x and ||^||r are 
non-zero outside S„ the calculations make sense point- 
wise at any point (r, i?) ^ E; all derivatives can be taken 
in the usual sense. If we can find such factors^S p and x 
in Ro with the further properties that neither x nor ||^||r 
vanish outside E„, and that equations (IMl) and (IMl) are 
satisfied away from E, then following the computation 
backwards we obtain (IMl) outside E. Then H'^x ~ ElpX 
belongs to C^{X), and we have shown that it is zero out¬ 
side E, and so simply zero since E has zero volume (is a 
set of zero measure). Thus Tpx = is an eigenfunction 
of H'. 

Let us now try a global treatment of the computation. 
In (IMl) . we used the Leibniz rule for derivatives. We do 
not know if it is valid here since it is possible that one 
factor contains a singular part (a non £^(A’)-part) which 
is compensated in another term. We would not be able to 
separate the terms but this is precisely what we must do 


in (IMli-dlZl). Assume that (IMl) is valid with each term in 
C^{X), possibly after restricting (^, x) to a smaller subset 
H' of Ho, (1301) . Now, we face the division problem in (1231) . 
Since p g H^{X), (H' — E)p is well-defined as a C'^{X)- 
function. So we should see an equality between C'^{X)- 
functions in (IMl) . Let us only consider x~^^”X ’ 
if one views as a distribution, there is the problem 
that the product of distributions is not generally defined, 
if one can identify x~^V"x with a distribution. Instead, 
we might view V^x ^nd as square integrable func¬ 
tions. Then V”x ’ is an integrable function in the 
R variable. But the multiplication by x~^ may destroy 
this integrability property. Another try could be to see 
X~^^”X ■ as a differential operator but, since x niay 
vanish, it would be a singular one. Again, the result of 
its action on p may be outside C'^{X). Anyway, we see 
that one already has difficulties even to give a meaning 
to (|M|). 

In Hunter’s formulation^, the first step of the compu¬ 
tation dUD - dMD is performed for a special choice of Xj 
given by (O . The full computation is followed in the pa¬ 
pers by Gross et al^—. They do not require the nuclear 
wavefunction to be square integrable and so the removal 
of the centre-of-mass motion is not performed. However 
it is perfectly possible, and convenient, to discuss their 
method in the framework set out here. The nuclear func¬ 
tion X is chosen as 

Y{R) = |vl/(r,i?)|2 dr)'^', (31) 

where S is an arbitrary real-valued function (cf ((3])) and 
T is a normalized solution of m- Then, p is defined 
by T/x and one derives equations for p and x as above. 
Recall that they interpret ^ as a i?-dependent electronic 
wavefunction. The meaning of (1221) already requires some 
information on the regularity of the factors. One can 
adapt the arguments i n^^i^^~ — to show that the function 

R J |'I'(r,i?)p dr 

is actually real analytic outside E„. Since it does not 
vanish there, its square root is also real analytic and so 
isxin (EU, if S is chosen real analytic. Away from E, the 
function p defined by 'L/x is therefore real analytic, by 
virtue of (EOI)- Thus we can follow our ‘local’ treatment 
performed above but, this time, we get a stronger result. 
The problem of factorisation of an eigenfunction of H' is 
equivalent to finding a solution {<p,x) of (l26l) and ((Ml) 
away from E such that ^x belongs to H^{X) and such 
that neither x nor ||^||r vanish outside E„. 

The difficulties described above in the global approach 
also appear in the work of Gross et a/.— as we now show. 
Of course, we have ^ but, since the integral in (IHTl) 
may vanish we have to be more precise in the definition of 
p. In view of (1211) and repeating the argument at the end 
of Section m we see that the integral in (1511) can vanish 
only in E„; thus p is well-defined outside E„. This would 
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be sufficient to define Tp everywhere as 'I' /x if it were in 
£?{X'). But the latter property is not certain since we do 
not know the behaviour of tk/y at the collisions, nor do 
we know if it is small enough at infinity for Tp to be square 
integrable. The same remarks apply to its derivatives of 
course. 

If one replaces the Coulomb interaction by a real an¬ 
alytic potential, then one can show by elliptic regularity 
that is real analytic everywhere and the above discus¬ 
sion is valid with an empty set S. So the computation 
outside S is actually the global one. Nevertheless, one 
would still have to characterize the behaviour at infinity 
of the product since it must belong to if it is 

to be an eigenfunction of H'. 


IV. VARIATIONAL METHOD 


A classical way to find solutions of a partial differential 
equation is to introduce an appropriate functional such 
that its critical points are precisely the solutions of the 
given equation. Then one tries to find local extrema of 
the functional. Cederbaum proposed to follow this strat¬ 
egy and introduced a functional on functions Tp(r, R) and 
x{R) and on two real parameters, having in mind that 
the product Ipx for a critical point should give a nor¬ 
malized eigenfunction of ht. Cederbaum’s arguments do 
not prove the existence of such factorised eigenfunctions 
since he does not prove the existence of critical points. 
There are important difficulties that make the search for 
critical points a delicate matter as we now describe. In 
the following we investigate the variational approach to 
the factorization of eigenfunctions of the internal Hamil¬ 
tonian H'. 

Let E be the set of the eigenvalues of H'. Recall that 
is the space of wavefunctions (with centre-of- 
mass removed) of Ai„ nuclei and N,, electrons such that 
'k and its distributional derivatives up to second order 
are all square integrable. We recall (cf (l30l) l that Ro 
is the set of couples, (^, x) where Tp g 'H^{X), and x £ 
•^2(]^3Af„-3) jg ^ nuclear wavefunction, and such that the 
product Tp{r,R)x{R) belongs to 'H^(A’). Consider the 
functional r : Hq x — >■ M specified by^ 

((^,x); (A,/r)) +\{1- ll^ll^) 

+A^(l-Ilxllfl)- (32) 


Let Hi he a subset of Ho such that r is differentiable 
on Hi X R^, and (^,x) G Hi, (A,/r) g R^. Then the 
following computations are valid 


— = 2x*(H'-A)(/3x, ^ = 2x(((<3|(H'-A)(/9x)r-I^) ■ 

((</?, x); (A,/r)) is a critical point precisely when both 
terms in (1331) are zero and 


ll<PXll = 1, llxlli? = 


(34) 


In such a case, 

fi|xl^ = X*(^I(H' - X)Vx)r = (^|x*(H' - A)^)^ = 0. 

Since x is not identically zero, p, = 0. 

Let TZy = {(^1 ^); x(-R) 7^ 0}- Denote by the com¬ 
plement, that is the set of zeroes of x- By the second 
equation in cannot be the whole space X. In¬ 

side Z^, the product Tpx is zero so that (H' — A)^ X = 0 
there. On the region — A)^ = 0, since we con¬ 

sider a critical point. If we assume that the boundary of 
Z^ (or Hy) lisis zero volume we can show that Tpx is an 
eigenfunction of H' as follows. We know that it belongs 
to H^{X). Thus (H' —A)^ is well-defined and belongs to 
C^{X). The latter is zero on TZ^ on Z^- Given that 
the boundary of Z^ has a zero volume, (H' — A)^x = 0 
holds true in C?[X') and, thanks to the second equation in 
(El, ipX is an eigenfunction of H' and A g £. Using again 
elliptic regularity, Tpx must be real analytic away from E. 
As already pointed out, this implies that Z^ C E„, which 
has a zero volume. We thus have shown that, if we have a 
critical point ((^, x); (A,0)) of r such that the boundary 
of Zy has zero volume, then Ipx is an eigenfunction of H' 
with energy A. 

This variational method is related to the system dMD 
and (1^ . Let ((^, x); (A, 0)) be a critical point of r. We 
assume further that Z^ C Eji and that x is smooth out¬ 
side E„. We know that (H' — A)^ = 0 and rewrite this 
away from E as 

0 = ^(T„x)+xHeZ^-A^ X- (35) 

Defining Eg/ by A - (x)"nTnX) away from E„, we get 
(1^51) outside E. Multiplying (1351) by Ip*, integrating over 
the electronic variables, and dividing by x, we recover 
(l29l) away from E. Using again the definition of we 
arrive at dH, outside E. 

An obvious difficulty for the present variational 
method is the determination of the space Hi, where the 
critical points of r live. As discussed in ApDendix IVIII Cl 
we think that Hi is strictly included in the space Ho, 
which is already not easy to describe. Another difficulty 
is related to our assumption that the boundary of the 
set Zy has zero volume. Without this assumption, we do 
not see how to justify that, for a critical point, the corre¬ 
sponding product Tpx is an eigenfunction of H'. We can 
only hope that an appropriate study of the critical points 
or the solutions of the system (|26|) and (l29l) justifies this 
assumption. 


V. A SPECIAL FACTORIZATION 

In the paper by Gross et al^ the factorisation m is 
chosen as the starting point; we saw in Section IIIII that 
it could produce singularities similar to those in Ceder¬ 
baum’s approach. Here we present another factorization 
based on an important result by Agmon^iiS^^ that does 
not have this drawback. 


I. 
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Let 'I' be an eigenfunction of H' with energy E which 
is isolated in the spectrum (this energy condition is satis¬ 
fied for relevant situations in Chemistry). Then, one can 
deduce from Theorem 4.13 in Agmon^i that there exists 
c > 0 such that the function exp(c|(7’, R) belongs 

to Here |(r, i?)| denotes the norm of the vector 

(r, A), that is, {\r\^ + \R\Y^^- Setting (A) = {1+\R\‘^)^/^, 
we choose the nuclear factor in the form 

x{R) = aexp(—c'(i?)) with a > 0 and c > c'V2 > 0. 

(36) 

Then x G 77^(1^^^^" ^^), x is smooth, real valued and 
positive everywhere. Choosing a appropriately, we can 
ensure that ||x||fl: = 1. Now we define ip = d'/y. Since 
we have pointwise (R) < V^\R\ < ^/2\{r,R)\ then 

exp(c'(i?))|«'(r, A)| < exp(c|(r, A)|)|^'(r, i?)| 

and ip belongs to C?{X'). If H is a bounded subset of A, 
then ip actually belongs to But we do not know if 

ip G 77^(A). In other words, the (distributional) deriva¬ 
tives of ip up to second order are locally square integrable 
but we do not have enough control of their behaviour at 
infinity to ensure that they are globally square integrable. 

These properties are however sufficient to allow us to 
repeat the arguments of Section Hill on any bounded sub¬ 
set n of X, starting from 41 = ipx, for then all terms are 
in £^(12). The nuclear wavefunction x is smooth, and so 
the explicit function Eei (i?) given by (1^ is also smooth, 
and bounded, and (l26l) even makes sense globally in X, 
in the sense of distributions, and takes the form 


the computation of Section imi to get (I37p . We can also 
look for a local extremum of r' at {ip,X), since it must 
be a critical point. Thanks to the fact that x is smooth 
and non-vanishing, we avoid the difficulties encountered 
in Sections mi and El Indeed, we may forget about the 
collision set E (although it reappears in the regularity 
properties of the solution of (1371) 1. and the computations 
are (almost) elementary. In the variational approach 
based on r', (l38|) . we avoid completely the obscure spaces 
T-Lq and Hi and the troublesome assumption about the 
boundary of Z^- 

Note that, in the above description, we have a rela¬ 
tively free parameter, namely c', that occurs in the defi¬ 
nition of X in (IM)l . We could try to let r' depend on c' 
as well in the variational approach. We do not know if 
this option facilitates the search for critical points. 

VI. NORMALIZATION 

In this section, we discuss the normalization used in 
Gross et alA and in Cederbaumi. In particular, we rec¬ 
tify a modification proposed in Cederbaum’s Erratuni^. 
In the framework of Section IIIIl let us assume that 
'I'(r, i?) = Tp{r,R)x{R)- Cederbaum’s normalization re¬ 
quires, for all R, that 

ll^llr = J MeR)\^ dr = I. (40) 

This implies that, for all R, 


( 


H' + 


ri^c'R 

MR) 



-3iV„(i?)2))(^. (37) 


Note also that, this time, the multiplication by xiR) 
preserves the space M{X). So, if we can find a nonzero 
solution ip of (EZD, such that ip € 77^(A), we can reverse 
the computation in Section m to show that ipx is an 
eigenfunction of H' with energy E. 

We can also modify the variational method presented 
in SectionElin the following way. Setting x(7?) as in (1551) 
above, we consider the functional r' defined on 77^(A) xR 
by 


{ip, A) !-)■ ((/3x|H'(<^x)) + A(I - \\ipx\?) ■ (38) 

In contrast to the functional r, r' can be differentiated 
everywhere and 

— = 2x*(H'-A)<px, -Tw = l-||£Xir (39) 

dip d\ 

Now, since x does not vanish, at a critical point (</5, A), 
the product (H' — X)ipx = 0. Since (H' — X)ipx G £^(A) 
this shows that ipx is a normalized eigenfunction of H' 
with energy A, and that X € £. As above, we can redo 


J I^(a^)P dr = \xiR)\^ ■ 

Thus (1311) is satisfied and the factorisation is the same 
as in Gross et alA. Conversely, if we start with the lat¬ 
ter, we have already pointed out that 'ip may be very 
irregular near the nuclear collisions. Fortunately, as seen 
at the end of Section imi X can only vanish at the nu¬ 
clear collisions, that is, in a small region (a set of mea¬ 
sure zero). Thus 'ip can be defined and we deduce that 
(l40l) holds true outside the nuclear collisions E„. Note 
further that 'ip ^ C?{X') so it cannot be interpreted as 
a molecular wavefunction; instead it should be viewed 
as a 7?-dependent electronic wavefunction as in the BO 
approximation. The normalization (I40p does not essen¬ 
tially change our discussion on global computations in 
Section imi Now assume, as in our discussion outside 
the collisions set E, that we have found a solution (^, x) 
of (1551) and (1551) away from E such that Tpx belongs to 
77^ (A) and such that neither x nor ||^||r vanishes out¬ 
side E„. Then 4^ = ipx but (itol) might be false. The 
factorisation 

'b(r,77) = • \m\r{R)x{R) (41) 

does satisfy the normalization condition but now the 
function iplr, R)/\\ip\\r{R) might be very irregular near 





S. Thus, it is not clear that the normalization (l40l) can 
be satisfied. 

Assuming that (^, x); j^st as above, solves the system 
(Hi and (1^51) outside S, we can derive the non-linear 
system that the factors in (IdTI) should satisfy, correcting 
in this way the corresponding computation in^. We define 
Lp = Tp/\\(p\\r and X = xll^llr (instead of tp = ^||^||r- 
and X = x/ll^l|r in^)- Away from E, ||(/7||r- = 1 and we 
know that (H' — E)T^ = 0, thus (H' — E)(px = 0, since 
V’X = WX- By the computation ( 1 ^ - ( 1 ^ . we arrive, still 
outside E, at 

Kelip = Eel{R)ip (42) 

where 

Eei{R) = -X-'(T„x)+i?andK,, = H'-^V"x-V", 

and 

TnX = [E - {p\\^eip)r)x- (43) 

VII. DISCUSSION 

In the previous sections, we have reviewed two schemes 
for the factorization of eigenfunctions of the molecular 
Hamiltonian, that were described in the papers^^— . We 
have seen that many ambiguities appear not only in the 
computations but also in the meaning of the results. 
Based on mathematical results (well-known in the mathe¬ 
matical physics community), we have extracted the main 
ideas of these methods and implemented them in a co¬ 
herent framework, giving in this way a precise meaning 
to the statements and partially justifying the computa¬ 
tions. We also have provided in (|V] a new factorisation 
by using the methods in a different way. 

The main results can be summed up as follows, start¬ 
ing with the first method. If an eigenfunction 4' of ( 1171 ) 
can be factored into Tp{r, R)x{R) with sufficiently regu¬ 
lar factors, then the latter must satisfy the non-linear 
system of equations (HH) and (1211), outside the set of col¬ 
lisions E. If one prescribes x to be a marginal of d' (see 
(ED)), one is led to the same conclusion, the factors be¬ 
ing automatically regular enough. Conversely, if one has 
a solution {Tp, x) of the non-linear system away from the 
collisions that satisfies some further conditions, the prod¬ 
uct ^(r, R)x{R) is an eigenfunction. If one chooses at the 
outset X as in equation (|36|) one essentially gets the pre¬ 
vious results without caring about the collisions. 

The second method is of variational nature. Having 
in mind to factorize an eigenfunction 4* as ^(r, i?)x(i?) 
one introduces a functional, acting on functions Tp and 
X, that is defined on an appropriate but quite compli¬ 
cated space. The functional is chosen such that its critical 
points produce an eigenfunction d)* in the desired prod¬ 
uct form provided they satisfy some regularity condition. 
Furthermore, one can relate this approach to the pre¬ 
vious non-linear system outside the collisions. In order 


to progress one must better understand the functional t 
and the complicated sets Ho and Hi. In Appendix I VHI Cl 
we show that "Hi is not empty and is probably strictly 
included in 'Hq. A demonstration that the functional t 
actually has critical points is a delicate matter. Further¬ 
more we need an unpleasant assumption on the zero set 
of X to construct an eigenfunction from a critical point. 
When the factor x is chosen as ED one can also follow 
a similar variational method on a quite natural, simple 
space. 

For both approaches the results are quite limited. In 
the factorisation with x(^) given by a marginal, the other 
factor Tp is smooth away from the collisions but we can¬ 
not, as yet, exclude an irregular behaviour near the colli¬ 
sions, which would not be compatible with an interpreta¬ 
tion of their product as a wavefunction. Starting from a 
solution of the non-linear system, we do not know if the 
conditions required to prove that the product is indeed 
an eigenfunction are satisfied. In the variational method, 
we do not know if the functional has critical points; we 
even have difficulties to describe the space where we have 
to look for them. If we have such a critical point, we still 
need to check further properties to ensure that the prod¬ 
uct of factors is indeed an eigenfunction. 

The situation is a bit better when one requires x to 
be specified by (1361) . In this case, the non-linear system 
reduces to a linear equation for the other factor Tp (see 
( 1571 ) ') and we just have to find a solution of this equation 
in a natural space (i.e. 'H^(A’)). In this framework, the 
variational method is also easier to work with. We know 
that the functional has critical points but we do not know 
how to compute them. 

In the previous sections, we have seen that the set of 
collisions E plays an important role. This is due to the 
fact that the Coulomb interactions have singularities pre¬ 
cisely on this set. These collisions are responsible for 
most of the difficulties we encountered. If we regularize 
each Coulomb singularity, that is, make the replacement 
of each in the operator H' by a real analytic func¬ 
tion, then an eigenfunction 4^ is everywhere analytic. If 
we assume that 4'(r, i?) = x{R)Tp{r,R) or if we write 
this factorisation with x dehned by ED, we can show as 
above that x never vanishes. All the previous difficulties 
related to the behaviour of Tp near E and those of x(^) 
near E„ disappear after the regularization. Therefore the 
regularized model presented in Gross et alA and the ex¬ 
actly solvable one chosen by Cederbaunt^ are not capable 
of giving insight into the actual molecular case, because 
an essential ingredient is lacking from the proposed mod¬ 
els. 

Given a factorization of a particular eigenfunction 
4'(r, i?)i = f{R)i(j){r,R)i, a ‘non-adiabatic energy sur¬ 
face’ for the nuclei can be defined formally by integrat¬ 
ing out the electronic variables in the expectation value 
of the internal Hamiltonian in the state (fii, (cf (l28l) l 

U(i?), = J ^{r,R)*H\r,Rmr,R),dr (44) 
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where H' is the molecular Hamiltonian with the centre- 
of-mass contribution removed (see ©). Hunter showed 
that one can derive a ‘reduced Schrddinger equation’ for 
the nuclear function f{R)i that partners R)i in which 
U(i?)i appears as a potential energy contribution. 

However it is important to note that this is not a re¬ 
finement of the conventional Schrodinger equation for nu¬ 
clear motion on a PES, because here the energy E is fixed 
(it is the eigenvalue of the specihed eigenfunction 
One can require f{R)i to vanish at the collisions (i? = 0), 
but the behaviour at oo is not a priori assured. In the 
conventional adiabatic BO treatment E along with the 
nuclear wavefunctions are unknowns, and one finds in the 
well-known way that satisfying the boundary conditions 
at i? = 0, i? = oo to assure square integrability is only 
possible for certain discrete values of E, the molecular 
vibration-rotation levels associated with the PES. 

The pseudo-potential U defined in (HI)) as introduced 
by Hunter—, and studied computationally by Czub and 
Wolniewica22, is only defined in a purely formal sense 
since, as we have seen fSection IHIL the function f may 
be so irregular that the application of H' to it could be 
ill-defined. It would be sufficient however that, for fixed 
i?, (f) lies in the Sobolev space Ti,^ in the r variables, but 
we do not see what would guarantee such a property a 
priori. Our analysis shows that, outside the collisions, 
(j) is analytic though that is not enough to control the 
behaviour at large r. However, making use again of ar¬ 
guments ii>i^i2i“— , one might hope to show that, for fixed 
R away from the nuclear collisions, (j) does have the regu¬ 
larity R? in the r variables, so giving a meaning to (Hil) . 
and also that U is analytic in this region. 

Hunter thought it unlikely that a nuclear function 
f{R)i would have zeroes in view of his interpretation 
of it as a marginal probability amplitude function for 
the nuclear coordinates, but originally based his claim 
on an analysis of the Schrodinger equation for coupled 
harmonic oscillators^^. That problem is no real guide to 
the properties of the solution of the Schrddinger equa¬ 
tion for the Coulomb Hamiltonian not least because the 
oscillator Hamiltonian is separable, has purely discrete 
spectrum and is well-behaved at collisions. 

Nevertheless, in the chemical physics literature on fac¬ 
torization of molecular wavefunctions it has been argued 
in more general terms that a nuclear wavefunction de¬ 
fined as a marginal probability amplitude for an exact 
eigenfunction, dt, as in (O or dSU, is necessarily nodeless. 
The argument rests on the statement that there exists 
a set of orthonormal functions {tT„(r, i?)} that is ‘com¬ 
plete in the adiabatic electronic space’; for definiteness, 
assume the {cr„} are the eigenfunctions of the clamped- 
nuclei Hamiltonian, 

H' - T„ ^ 

in which the nuclear positions {R} are treated as classical 
parameters, 

Hcn(l^ ■ R)^n(^‘^: R') — ^n{R')^n(^‘^: R') 


Then one writes an exact eigenfunction (of H') as in the 
Born-Huang theory^S 

R) = Xn{R)(Tnir, R) (45) 

n 

and the nodeless property of the {xn} follows^d^i^l. 

By contrast we are unable to exclude the possibility 
of nodes in the x{R) functions for R values associated 
with the collisions. Why the difference ? The essential 
point is that an expansion such as (145)1 relies on the set 
of eigenfunctions {fTn} providing a resolution of the iden¬ 
tity, and this is only valid in the case of a purely discrete 
spectrum; in other words, {cr„} must be true eigenfunc¬ 
tions associated with (discrete) eigenvalues. When an 
operator also (or only) has a continuous portion of spec¬ 
trum the matrix notion of diagonalization providing a 
complete set of states breaks down, and the resolution of 
the identity must instead be developed from the spectral 
theorem and the idea of spectral projection. One of us 
has explored this idea— in detail in the context of the BO 
approximation; we refer to that discussion which shows 
that an exact representation of 41 is much more compli¬ 
cated than (HSl) . and does not lead us to such a definite 
conclusion about the nodal properties of %. 

So far, nothing has been said about spin statistics. 
Consider a collision involving two identical nuclei (1,2). 
Under a permutation P 12 an exact eigenfunction Ik will 
either be symmetric (boson statistics) or antisymmetric 
(fermion statistics). In the later case the eigenfunction 
vanishes at the collision R 12 = 0, and a % factor calcu¬ 
lated according to ([3]) will also vanish. Direct examina¬ 
tion of the Schrodinger differential equation in the vicin¬ 
ity of Ri 2 = 0 shows that the spatial part of the wave- 
function may vanish in any cas o^^d^ . 

Based on the present knowledge, we have the follow¬ 
ing impression of these methods. The system (l26l) and 
(1291) is non-linear and has a priori singularities. This 
is already a difficult problem, but here we have the un¬ 
usual situation in the case of (1351) where the singularities 
of the equation depend on an unknown function. Con¬ 
cerning the variational method, the space on which we 
can apply it is difficult to describe. This comes precisely 
from the fact that the functional contains products of 
the variables, which are a priori less regular functions. 
When one requires that the factor y or y is given by an 
appropriate exponential function, as in Section El the 
situation is a bit better but we do not see a real improve¬ 
ment in (1371) compared to the original equation (ITTl) . We 
do not see a natural physical interpretation for the fac¬ 
tors in this setting and judge the factorisation artificial. 
Eor these reasons, we are not convinced of the efficiency 
of the methods to produce eigenfunctions. It would seem 
that the equation systems suggested in these proposals 
are so difficult to handle that a direct approach to con¬ 
structing eigenfunctions treating the electrons and nuclei 
on the same footing might be no more challenging; after 
all, the equations for the ‘electronic’ factor (p (or Ip) still 
contain the full internal molecular Hamiltonian H' and 
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all the electronic and nuclear variables of the problem. VIII. APPENDIX 

In this Appendix we explain some notions and results 
used in the main text. 


A. DifFerentiation in the distributional sense and products 

Let d be an integer and / be a locally integrable func¬ 
tion on R'^. This means that, for any bounded subset 
n oiW^, f is integrable on fl. We denote by the 

space of smooth (complex-valued) functions on R*^ with 
bounded support. One can identify the function / with 
the distribution Tf : T>(R‘^) —>■ C (the set of complex 
numbers) defined by 

Tf{g) = [ f{x)g{x) dx 

The distributional derivative of / w.r.t. xi is the corre¬ 
sponding derivative oi Tf, which is the new distribution 
dx^Tf : 'D{W^) —>■ C defined by 


9x^Tf{g) = - f{x)dx,g{x) da;. 
aRrf 

Note that, for g G 2?(R‘^), dx^g G T>(R^). If h is a smooth 
function on R"* , then hg G V{R'^), if g G If T 

is a distribution on R'^, that is a continuous (in an ap¬ 
propriate sense) linear map from T>(R'^) to C, one defines 
the product of T by the smooth function h as the new 
distribution given by {hT){g) = T(hg). In particular, 
hTf = Thf. 

Let us take examples that are relevant for the main 
text. For simplicity, we assume Ne = Nn = 1. Let 
x(i?) G L^(R^) and ip{r,R) G L^(R®). The distributional 
derivative of x w.r.t. i?i is the linear map 

dniX ■ T>{R^) 5 g ^ - f x{R)dR^g{R) di?. 

Jr? 


We may multiply it by a smooth function. But the prod¬ 
uct by g}{r, ■) (with fixed r) is a priori undefined since 
T{x,-)g G T’(R^) could be false. Indeed <p{r,-) could be¬ 
long to the set of L^(R^)-functions that are not smooth. 
So, we have difficulty defining Lp(r,-)dR,^x- However, if 
Br^X = T/j, for a function /i G L^(R^), we may define 
the product ;p(r, •)9 _RiX the usual product (p{r,-)fi, 
which is in L^(R^). We could also differentiate this prod¬ 
uct in the distributional sense since one can identify a 
L^(R^)-function with a distribution. If we would do that, 
the Leibniz rule 

dR,{T{r,-)fi) = {dRMr,-))fi+T{rr){dRji) 

might be false. Both products on the r.h.s might be un¬ 
defined. This simple situation illustrates some difficulties 
mentioned in Section imi 
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B. Products of "H^-functions 


since (px = 0. Now, for j large enough, 


Here we construct functions tp € 7^2j-]^3(7Ve+7V„_i)^ 

X G 'H2(K3Ar„_i) ^ ■h 2(IR3(^=+^—i)). Let 

/ and g be smooth functions with bounded support in 
^3(Ne+N„-i) ]^3Ar„_i^ respectively, such that both 

are equal to 1 near 0. Let a and /3 be real numbers. We 
set 

p(r,Ii) = /(r,i?).(|r|2+|i?|2)“/2andx(i?) = gWW ■ 
Then p G 7^2^]g3(Afe+iVTi-i)^ and x G if 

a > 2 - ^(iVe + iV„_i) and /3 > 2 - ^^n-i ■ 

Provided that iVg > 3, one can choose a and /3 satisfying 
the above conditions and also a + /3 < —3A^„_i. The 
latter implies that px ^ £2(]R3(Afe+Afn-i)^ and thus px ^ 
-^2(]^3(Are+JVn-l)^^ 


C. Hi ^ Hoi 

To begin with, if p and x are smooth functions with 
bounded support in M3(Afe+Af„) respectively, 

then {p,x) G ^1- Thus Hi is not empty. By defini¬ 
tion, Hi G Hq but even so we expect that one can find 
{p, x) G Hq at which the functional r is not continuous 
(and thus not differentiable). To motivate this guess, we 
shall prove it for a simpler functional related to r (Sec¬ 
tion nYi). Consider the map tq : Hq —?► M given by 

which is a part of the functional r. Take Nn > 3. Let 
P = 2- 3Nn/2 + 1/8 > 2- 3iV„/2. Let ^ = 3A^„/2 -f 
/3 -I- 1/8 = 2 -I- 1/4. In particular, 2 + 5 < 3iV„/2. Let 
/ be a non-zero, smooth function on R^JVe with bounded 
support. Let g and h be two smooth functions on 
with bounded disjoint supports such that g = 1 near zero. 
Since the supports are disjoint, gh is identically zero. Let 

+{r,R) = f{r) ■ h{R) and x(i?) = giR) ■ \Rf ■ 

Thanks to /3 > 2—3Nn/2, {p, x) belongs to Hq; obviously 
PX = 0 identically. For all integer /, let gj{R) = j^g{jR)- 
Notice that, for j large enough, the support of gj is in¬ 
cluded in the region about 0 where g = 1. Let 


Pj{r,R) = f{r)-gj{R). 

Since 2+5 < 3N„/2, we see that gj goes to 0 in 

and thus pj goes to 0 in as j goes to 

infinity. 

If To were continuous at (i^, x) then the difference 
TQ(p+pj,x)~To{p, x) should go to 0, as j goes to infinity. 
This difference is given by 

TQ{p + Pj,x)-To{+,X) = W+jXf , 


= ll/ll^ / \gj{Rr-\R\^^dR 

jRSNn 

^^.2.-3M||^||2. f |^(,)|2.|g|2/5dg 
Jr 3 N„ 

and the r.h.s. blows up as j goes to infinity, since 25 — 
3iV„ —2/3 = 1/4 > 0. This yields a contradiction, showing 
that To is not continuous at {p,x)- 


D. (/> not necessarily in H^ 


The point we want to demonstrate here is simply that 
one must make a check to ensure that a factor (j)(r;R), 
defined by (O, belongs to the Sobolev space H^{r). This 
is an essential property if p is to be interpreted as a wave- 
function for the electrons for fixed R. Here we propose 
two model examples that show that 4'// is not automat¬ 
ically in H^ (in the r variables) when 4* belongs to H^ 
(in all variables, r, R) and / is defined as a marginal as 
in ([3]). In the first model 4'// does not belong to H^ 
because of irregular behaviour near i? = 0; in the sec¬ 
ond the irregular behaviour is located at oo (in the R 
variable). 

We choose a smooth, nonnegative function t of one 
real variable such that 


r(t) 


1 if |t| < i 
0 if |t| > 1 


For simplicity we choose three dimensional variables r 
and R with r = |r|,i? = |R|. The configuration space 
H = R®. Let 


4'(r,R) = R” t{R) T{R^"^{r - R)) 

for some integers n > 0,m < 0. One can specify n,m 
such that this ‘wavefunction’ 4* belongs to the Sobolev 
space 'H^(H) (essentially |n| has to be large compared to 
|m|). 

Following the prescription of Hunter— we then have, 
after integration over the angles and some simplification, 
for some c > 0, 


fiR) = cR”+i-™T(i?) • 


1 


[ T{tf {l+tR-^-^- 

Jm. 



Note that the vanishing of / at i? = 0 is stronger than 
that of 4'. Now we define 

f 0 if t{R) = 0 else 
(j>[r; R) = < 

i ~7wr 

In particular, we do not define (j) on the zero volume 
region {i? = 0}. Using the same changes of variables 
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as above, we get for some c' > 0, in the region where 
t{R) ^ 0, 

(j){r;R) = c'R^-^T{R^^{r - R))- 

JR / 

Explicit calculation shows that the ‘electronic’ function 
4> is square integrable, but that its first derivative (in 
the r-variable) is not for \m\ large enough, so that it no 
longer belongs to . Nevertheless there are n values 
such that d' is in ( n > —7m/2 — 2). The reason 

for this behaviour comes from the fact that / vanishes 
more strongly at i? = 0 compared to d'. 

We can use this idea again to translate the irregular 
behaviour at i? = 0 to i? = oo . Taking 

«'(r,R) = i?” (l-r(i?)) T(i?2”^(r - R)) 

but now for large —n and positive m, we can check that 
dt belongs to the Sobolev space (essentially \n\ has 

again to be large compared to |m|). Again we can adjust 
m such that the r-gradient of cj) is not square integrable. 

E. Local Sobolev space on R3{Afe+Ar„) 

In the text, we defined the Sobolev space 
•^2(]^3(Afe+Afn))^ Similarly, we can define 7^^(fl), for any 
bounded open subset 17 of The correspond¬ 
ing local Sobolev space, denoted by is 

the space of those functions / that belong to 77^(17), 
for all bounded open subsets 17 of A 

function that has usual continuous derivatives up 
to second order always belongs to the local Sobolev 
space but the integral of its modulus 

square can be infinite (just think about the constant 
function equal to one). Therefore the local and global 

■H2(R3(Afe+iV„)) 

-spaces are different, the latter being 
included in the former. 


F. H has no eigenvalue 

Let us first give a ‘physical proof’. Since the full molec¬ 
ular system is considered as being isolated in the uni¬ 
verse, its mass centre is freely moving. Thus it cannot be 
in a bound state. 

Now we turn to a mathematical proof. We start with 
d?]) (assuming for simplicity that h? /2 Mt = 1) and recall 
that the internal Hamiltonian H' is R-independent. Now 
we view this formula in the Fourier space of R € 
Denoting by Hi this representation of H, we get Hi = 
H' -I- M/, where M/ is the multiplication operator by the 
function f{^) = (^ being the Fourier variable asso¬ 

ciated to R). This can be rewritten with the help of a 
direct integral (cfi^, p. 279-287) as 

Hi = + H') de 


According to^, p. 284, E is an eigenvalue of Hi (or H) 
if and only if the 3-dimensional volume of Ce = {C S 
Ais an eigenvalue of|^p + H'} is positive. Note that 
E is an eigenvalue of + H' if and only if A — 
is an eigenvalue of H'. By the Mourre theory (see^), 
one can show that the set of eigenvalues of H' is at most 
countable. This implies that Ce is the union of an at 
most countable set of spheres. Thus Ce has zero volume 
and H has no eigenvalue. 


G. Elliptic regularity 

We propose here a short and intuitive introduction to 
elliptic regularity. For more details, we refer to^i^. 

We consider the differential equation = / on 
where the function / is given and u is the unknown func¬ 
tion. Our goal is, knowing the regularity properties of /, 
to obtain those of any solution u. In one dimension, the 
problem is easy, since the second derivative of u is exactly 
/. In the present case, it could happen that the deriva¬ 
tives d'^u and dyU have singularities that cancel when the 
sum is performed. 

A good way to study the regularity of u is to consider 
the Fourier transform ft of tt. Indeed, regularity proper¬ 
ties of u are encoded in the decay properties at infinity of 
u. This can be seen from the following (formal) identity: 

Vu(x) = (27r)“^ f ipu{p) dp. 

7r2 

If u has a fast enough ‘decay’ in \p\ at infinity, it makes 
the above integral absolutely convergent and u' is nice. 

After Fourier transformation, the equation becomes 
\p?u{p) = f{p). For IpI > 1, we get u{p) = \p\~‘^f{p)- So 
the ‘decay’ of u is better than that of / and u is more 
regular than /. This is called elliptic regularity. 

Now we can also apply this method to the equation 
V'^u+Vu = / on where H is a function. For instance, 
if V and / are smooth, we can write u{p) = \p\~'^g{p), 
where g = f — Vu. Given some regularity for u, we 
can improve it by the previous formula. The improved 
regularity can be plugged into the formula again to get 
a better regularity, and so on. If H or / has a limited 
regularity, so does u. 

Sometimes one is forced to view the equation -I- 
Vu = f in the distributional sense. In this case, one 
can still follow the above argument and, when / and V 
are smooth, so is also u, and u satisfies the equation 
+ Vu = f in the usual sense. 

In the main text, we used the elliptic regularity to 
get the real analyticity of the solution of an equation. 
Roughly speaking, a function u on is real analytic if, 
near any point (xq, yo), it can be written as a polynomial 
of infinite degree (i.e. a series in powers {x—zoYip—yoy). 
A real analytic function is always smooth but the func¬ 
tion (a:, y) !->■ xe~^l^, if ce > 0, and (x,y) >->■ 0, if a: < 0, 
is actually smooth but not analytic. 
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If V and / are analytic, then so is any solution of 
V^u + Vu = f. The argument to see thisii is more 
involved than the one above. In the Coulomb case, the 
potential is real analytic away from the set S of collisions. 
So the above arguments apply away from S. If you think 
about the hydrogen atom and set r = x'^ + y'^ + z'^, the 

function {x, y, z) i->- is a solution of + M(l/r — 

1/4) = 0. The potential and the solution are real analytic 
away from 0 (and the latter has to be since the equation 
is elliptic) but they both are not even smooth at 0. 
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